Superfluid Suppression in rf-Wave Superconductors due to Disordered Magnetism 
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The influence of static magnetic correlations on the temperature-dependent superfluid density 
paiT) is calculated for d-wave superconductors. In self-consistent calculations, itinerant holes form 
incommensurate spin density waves (SDW) which coexist with superconductivity. In the clean limit, 
the density of states is gapped, and p s (T -C T c ) is exponentially activated. In inhomogeneously- 
doped cases, the SDW are disordered and both the density of states and p s {T) obtain forms indistin- 
guishable from those in dirty but pure d-wave superconductors, in accordance with experiments. We 
conclude that the observed collapse of p s at x ~ 0.35 in underdoped YBa2Cu306+:r may plausibly 
be attributed to the coexistence of SDW and superconductivity. 
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I. INTRODUCTION 



High temperature superconductors (HTS) are an ideal 
class of materials with which to study electronic corre- 
lations in superconductivity because the correlations can 
be tuned from weak to strong via chemical doping. A 
consequence of strong-correlation physics is that the BCS 
theory of conventional superconductors fails to describe 
HTS. Uemura demonstrated that, unlike in BCS theory 
where the superfluid density p s and critical temperature 
T c are independent, HTS exhibit an approximate scaling 
p s oc Tc 1 . The physical origin of the Uemura relation 
is not known conclusively but is consistent with strong- 
correlation models 2 * 3 -^ in which the quasiparticle spectral 
weight, and consequently p s , are proportional to the hole 
concentration p, where p is measured relative to the Mott 
insulator phase. However, these models fail to explain 
both the collapse of superconductivity at a nonzero dop- 
ing p c ss 0.05 and the breakdown of the Uemura relation 
near p Cl shown by recent experiments in YBa2Cu306+ x 
(YBCO)££. 

A number of authors^ have suggested that T c is 
governed by phase fluctuations, possibly in combina- 
tion with quasiparticle excitations^, and in particular 
that the rapid collapse of p s and T c at p c can be thus 
explained 10 .. In this work, we examine a completely dif- 
ferent mechanism: the suppression of superfluidity by 
the formation of static magnetic moments. This is mo- 
tivated by a substantial body of experimental evidence 
for the presence of quasistatic magnetism in underdoped 
HTS-iii^ii^ii&iLi&i 9 . There have been previous sug- 
gestions that some form of competing order is impor- 
tant in the underdoped HTS ) 20 ' 21 ' 22 ' 23 and in particu- 
lar there have been numerous studies of the competi- 
tion between <i-wave superconductivity and commensu- 
rate antiferromagnetismi 24 i 25 i 26 i 27 i 28 These are generally 
difficult to reconcile with superfluid density measure- 
ments largely because the competing order introduces 
an identifiable energy scale. Calculations show that this 
energy scale appears in the temperature dependence of 
the superfluid densit y 29 ! 30 ' 31 but such an energy scale 
has not been observed experimentally 31 . Indeed, recent 
microwave conductivity measurements 6 of the superfluid 



depletion, Sp s (T) = p s (0) — p s (T), in high quality single 
crystals of YBCO find 5p s (T) oc T with a crossover to 
5p s (T) oc T 2 when T < T c . The linear T-dependence 
is expected in a single-phase <i-wave superconductor and 
the low-T crossover to T 2 behavior has been attributed 
to residual impurity scattering. It is therefore not a pri- 
ori clear that the experimentally observed magnetic mo- 
ments have any significant effect on the electronic spec- 
trum. Here, we show that a phase of coexisting spin 
density wave (SDW) and <i-wave superconducting (dSC) 
order can, provided the SDW is disordered, have a spec- 
trum indistinguishable from that of a dirty dSC. We con- 
clude that the rapid collapse of superfluid density near 
p c could indeed be due to magnetism. 

Our approach is semi-phenomenological. We con- 
struct a mean-field model in which the model pa- 
rameters are assumed to have been dressed by elec- 
tron interactions. The approach is motivated by a 
variety of calculations ^i 3 - ' 24 ' 32 ' 33 ' 34 mostly for the t-J 
model, in which mean-field theories are developed for 
which the parameters are functions of p. In the sim- 
plest Gutzwiller approximation for the t-J model^ for 
example, the renormalized kinetic energy operator T 
is related to the bare kinetic energy operator To by 
T = To2p/(l + p). Other results are found in other 
approximations ; 2 ' 3 ' 24 ' 32 ' 33 but all show the same qualita- 
tive result that T is reduced as one approaches the Mott 
insulating phase. In dynamical mean-field theory calcu- 
lations, this derives from a self-energy which renormal- 
izes both the quasiparticle spectral weight and effective 
mass. 35 We remark that the effective interactions are also 
expected to depend on the doping, generally increasing 
as p is reduced. This is ignored in our calculations since 
it will have a quantitative but not qualitative effect on 
the outcome. The essential physics in the current discus- 
sion is that, near the magnetic phase boundary a small 
change in the mean-field parameters produces a much 
larger change in the magnetic state. The progression 
from pure dSC to pure magnetic order depends only on 
the general trend that the ratio of kinetic to interaction 
energies decreases as p decreases. 

In section [TT1 we introduce the model and describe the 
phase diagram. The most significant result of this section 
is that it is possible to have substantial incommensurate 
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magnetic moments coexist with the superconductivity, 
with very little suppression of the pair amplitude. This 
is in contrast with the more widely studied case of com- 
mensurate magnetic order, which suppresses supercon- 
ductivity rapidl y 24 i 25 ' 26 i 27 ' 28 In section IIII1 we calculate 
both the density of states and the superfluid density for 
the incommensurate phases. As mentioned previously, 
we find the surprising result that when the magnetic mo- 
ments are disordered, the spectrum is indistinguishable 
from that of a dirty d-wave superconductor. We conclude 
briefly in section llVl 



II. MODEL AND PHASE DIAGRAM 

The HTS consist of conducting two-dimensional CuC>2 
layers that are weakly coupled along the perpendicular 
direction. We model the lower Hubbard band of a single 
two-dimensional layer with an extended Hubbard model, 
treated at a mean-field level. Our numerical calculations 
have found that the results are only weakly dependent on 
the filling p but depend sensitively on the quasiparticle 
bandwidth and Fermi surface curvature. As discussed 
above, we assume that doping effects occur indirectly 
through a parameter w(p) which renormalizes the quasi- 
particle dispersion. For simplicity, and since the detailed 
relationship between w and p is not established, we will 
treat w as the independent parameter, and keep all other 
parameters fixed. 

Calculations are for an AT-site two-dimensional tight- 
binding lattice with periodic boundary conditions and 
lattice constant a = 1 , similar to one used previously to 
study the local density of states in underdoped HTS 36 . 
The Hamiltonian is 



E 

(i,3) 



. (1) 

where fcy are the hopping matrix elements of the tight- 
binding band, i and j are site-indices, a is the electron 
spin, and a = —a. We take ty = to for i = j and 
tij = t n , n = 1,...,3, for nth nearest-neighbor sites i 
and j. Taking {to, ...,t 3 } = {1.7, -1, 0.45, -0.1} gives 
the Fermi surface shown in Fig. [T] The local electronic 
density rii a = (hi a ), where hi a — c\ a Ci a , is determined 
self-consistently and the hole density is pi„ = \—rn a . The 
magnetic order parameter is then m; = (n^ — n^)/2, and 

the staggered moment is mf = {—l) Xi+Vi nii where = 
(xi, yi) is the coordinate of site i. The nearest neighbor 



pairing term A, 



-(fij), with/ij = (c 3i c lT -c jT c a )/2 



is also determined self-consistently and has pure dSC 
symmetry, Ay = ^(— \) y i~ Vi , in the nonmangetic phase. 

To allow for inhomogencous doping we add dopant- 
impurity and Coulomb interaction terms, also treated at 
the mean-field level: 



N, 



H, 



v j )h i n j -Z^y^V{T i --Ri)h i , (2) 
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FIG. 1: (Color online) Phase diagram as a function of kinetic 
energy renormalization. Open symbols are for the homoge- 
neously doped calculations, closed triangle and + symbol are 
for inhomogeneously doped calculations. Inset: Fermi surface 
and CSDW nesting vector. 
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with m — ~Yli(j n i<y-, Z the impurity charge, R,£ the loca- 
tions of the ATj impurities, and V(r) = e 2 e~ r / A /ea,Q a 
weakly-screened Coulomb interaction, with e 2 /eao = 1.5 
and A = 20ao- In the cuprate HTS, donor impurities 
typically sit a few A above the Cu02 layers, so we ran- 
domly choose values of = [xt,yt,d) with d = ao- 
The total hole doping is p — ZNj/N. We study a ho- 
mogeneously doped case with Nj = N and Z — p, and 
an inhomogeneously doped case at the same filling with 
Nj = N/4 and Z — Ap. The resulting impurity potential 
is smoother than one expects in many underdoped HTS 
but is reasonable for underdoped YBCO where approxi- 
mately 35% of chain oxygen sites are filled. It is assumed 
that strong-correlation renormalizations of V, J and U 
are included implicitly and remain constant over the nar- 
row doping range explored here; for a given p, we choose 
J and U such that w = 1 corresponds to a pure dSC 
phase close to the magnetic phase boundary. We then fol- 
low the magnetic phase diagram, Fig. [TJ as w is reduced. 
The results depend sensitively on the Fermi surface shape 
(ie. on t% and t^), but depend only weakly on p which is 
therefore chosen for computational convenience. We have 
studied the parameter sets (p, U, J) = (0.05, 3.4, 1.8) and 
(p,U,J) = (0.35,3.2,1.5). The two agree semiquanti- 
tatively where we have been able to compare; however, 
it is difficult to obtain converged solutions for p = 0.05 
when w is small, and we have chosen to present results 
for p = 0.35 where a full set of results is available. 

The calculations proceed as follows. The Hamiltonian, 
written H = c^Hc with = [eL . . . c]y»cij . . . cjyj.] and 
H a Hermitian matrix, is diagonalized numerically giving 
the unitary matrix U of eigenvectors and eigenvalues E n . 
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FIG. 2: (Color online) Typical self-consistent solutions. The 
staggered magnetization is shown for the CSDW phase at 
w = 0.85 (a,b) and the DSDW phase at w = 0.8 (c,d) with 
homogeneous (a,c) and and inhomogeneous (b,d) doping; The 
magnetization (e) and dSC gap (f) are shown for w = 0.78 
with inhomogeneous doping. Inhomogeneous results are for a 
single dopant configuration. 



The calculations 

n,E n <0 



mi 



— 1 — 2^ \^i+Nn\ 



(3) 
(4) 



n,E„<0 



\j = -(<^/2) 2J (U in U j+Nn + U i+Nn U jn ) (5) 

are iterated until the largest difference between succes- 
sive values of and n icr is less than 10~ 4 . At con- 
vergence, the total energy is typically varying in the 
tenth significant figure. Convergence is difficult to ob- 
tain: first, the charge density oscillates wildly in simple- 
iteration schemes because of the Coulomb interaction 
and, second, the magnetic moment configuration may fail 
to converge because the energy near self-consistency de- 
pends only weakly on it. A significant effort has been 
made to address both these issues. First, we have gen- 
erated our initial guess for the local charge density us- 
ing a self-consistent Thomas-Fermi calculation. In order 
not to bias the outcome of the calculation we have used 
8-16 randomly seeded initial magnetic moment config- 
urations for each parameter set, from which the lowest 



energy self-consistent solution is retained. Second, we 
have adopted a Thomas-Fermi-Pulay iteration scheme^ 
which controls the iteration instability in most cases. Re- 
sults shown here have all converged. 

The phase diagram is shown in Fig. [TJ For large w, 
there is a pure dSC phase. As w is reduced, there is a 
second order transition into a coexisting phase of dSC 
and checkerboard SDW (CSDW) order at w rj 0.95, fol- 
lowed by a first order transition into a phase of coexist- 
ing dSC and diagonal SDW (DSDW) order at id 0.81; 
both phases are illustrated in Fig. O Superconductiv- 
ity is destroyed at w ks 0.76. From Fig. [21 one sees 
that doping-induced inhomogeneity disorders the SDW 
and has a significant effect on the phase diagram in the 
dSC+DSDW phase: a typical solution for w — 0.78, 
shown in Fig. [2] (e,f), consists of an inhomogeneous mix- 
ture of pure SDW and dSC+SDW order, with supercon- 
ductivity preferentially forming in hole-rich regions. The 
spatially-averaged A varies considerably between dopant 
configurations, as seen in Fig. [TJ but within supercon- 
ducting domains, the local order parameter is consis- 
tently |Ajj | ~ 0.5; the destruction of superconductivity 
occurs inhomogeneously. 

We can understand the origin of magnetic order in the 
CSDW phase; the Fourier transform m(q) of rtii has a 
set of four peaks at (tt ± 8, tt ± 8) and the inset to Fig. [JJ 
shows that the vector (tt — 8, tt — 8) taken from the data 
in Fig.^a) connects nodal points on the Fermi surface^. 
The CSDW, therefore, nests portions of the Fermi sur- 
face where the pairing energy is small, and consequently 
minimizes competition between magnetic and dSC or- 
der. This explains why A is roughly constant through- 
out the dSC+CSDW phase (c.f. Fig. [TJ, and why the 
transition to the dSC+DSDW phase appears to occur 
only when the magnetic energy scale M ks U max(|mj|) 
is greater than A. It appears, then, as if CSDW ordering 
is stabilized by superconductivity. An analysis of m(q) 
for the DSDW, in contrast, does not reveal any nesting 
of high- symmetry points, and A collapses rapidly in the 
dSC+DSDW phase. 

These results arc in striking contrast to what one 
finds for the case of commensurate magnetic or- 
der. We show, in Fig. [3] the results of calcula- 
tions for the antiferromagnetic moment m z and the 
dSC order parameter A determined self-consistently 
in the homogeneous limit. The calculation proceeds 
as follows: Adopting a four-component notation 26 
A k = (ckf, cl k _p c k+QT , &_ 

-1/2. 



where c klJ = 
r,), one can write H = 
Ak^HkAk where the prime indicates a sum over 
k x > and 
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(6) 



with the band energy ek 
4^2 cos k x cos k v and Ak = 



= t + 2ii(cosfc x + cosky) 
^(cosk x — cosky) and M 
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FIG. 3: (Color online) Phase diagram for commensurate or- 
der. Here, U = 3.6 while other parameters are as given in the 
text. 



Um z . If the 4x4 matrix diagonalizing Hk is denoted 
Uk, then the self-consistent equations for A and m z are 



A = E ^2 (cos Mr - cosk y )Ui jk U 2j k, (7) 



Ni 



3 = 1 k, 

E jk <0 



x 4 

z = ^;E E [^ijk^sjk + u 2 jkUijk]. (8) 



j=l k, 

E jk <0 



Figure [3] shows that the dSC and antiferromagnetic order 
parameters are generally incompatible, with only a small 
coexistence region. More extensive studies of the phase 
diagram with commensurate order by Kyung2£ show that 
the size of the coexistence region depends on the model 
parameters, but that the antiferromagnetic and dSC or- 
der always suppress one another. By contrast, the CSDW 
order has very little effect on the dSC phase. 



III. DENSITY OF STATES AND SUPERFLUID 
DENSITY 



Because of the nodal nesting, the density of states 
(DOS) 

N 2n 

N ( u ) = EEH^I 2 ^ - S «) + \Ui+Nn\ 2 S{uJ + E n )}, 
i—l n—1 

(9) 

develops a gap of width 8 in the dSC+CSDW phase, 
as shown in Fig. Q] In contrast, the dSC+DSDW phase 
retains the characteristic d-w&ve DOS, N(uj) oc \oj\ at low 
u>, but develops a resonance at ujq which in many cases 
dominates the spectrum. In both cases, these qualitative 
differences from the pure dSC DOS are reflected in p s (T). 



FIG. 4: (Color online) Density of states. N(u>) is shown 
for homogeneous (top row) and inhomogeneous (bottom row) 
doping. Inhomogeneous results are for a single dopant config- 
uration, except for w = 0.78 which is averaged over 5 config- 
urations. The single-particle gap 8 and subgap resonance at 
uio are indicated. 



The superfluid density is related to the magnetic pene- 
tration depth A(T) measured in experiments by p s (T) = 
mc 2 /4ire 2 \ 2 (T). In linear-response theory, A~ 2 (T) = 
(Aire 2 /c 2 )(K^(T) - K^(T)) a ^ y with 



= Efc 1 ] mm /(^) 



(10) 



T^para 



r~ i r~ i f( E «i) - f(E n ) 



where 

pvr 1 ] 



/ ' L ' Q J ran L <P\ nm T? _ TP 



i,j p=0 
1 

E E U m i+pN Yl^h Uj+pN n 
i,j p=0 



where M and 7 are the inverse effective mass tensor and 
current vertex respectively, written in the basis of eigen- 
states of the Hamiltonian. On a tight-binding lattice, 



\M 



a/3 nj 



and [7 Q ]y = ic 



r^tij with 



a and P the unit vectors x or y and r.j 



The 



calculations are restricted to low T where we can use the 
T = values for A^ and ; the T-dependence of p s (T) 
is due to thermal pair breaking, as has been argued in 
Ref. [2|. We discuss this assumption below. 

We focus first on homogeneous doping, Fig.[5fa). The 
superfluid density at T = is a strong function of w, 
especially in the dSC+DSDW phase; p s (0) is reduced to 
10% of its initial value with only small changes in the 
coherence peak energy in the DOS. While this is con- 
sistent with experiments, the T-dependence of p s {T) is 
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FIG. 5: (Color online) Superfluid density. p 3 (T) is plotted for 
homogeneous (a) and inhomogeneous (b) doping. Different 
curves for w = 0.78 in (b) correspond to different randomly- 
generated impurity potentials. Inset: low-T behavior of p s (T) 
for homogeneous (solid) and inhomogeneous (dashed) doping. 
Results in (b) are configuration-averaged over 2-4 samples, 
except for w — 0.78 where results are shown for each config- 
uration. 



moment: however, the paramagnetic response at T = 0, 
i4T para (0), is a strong function of the magnetic moment. 
This is reminiscent of the response to disorder in d-wave 
superconductors where Cooper pair breaking by impurity 
scattering manifests as a nonzero JC para (0). In this case, 
however, the broken Cooper pairs are also apparent as a 
finite residual density of states at the Fermi level. A dis- 
order level sufficient to cause a 90% reduction in p s (0), 
as we have found here, would produce a residual density 
of states comparable to that of the normal state. The 
fact that such a residual density of states is not observed 
in our case illustrates that the SDW correlations are not 
simply breaking Cooper pairs. 

Rather, it is the fact that Cooper pairs in the magnetic 
phase do not have a well-defined charge-current which is 
responsible for the suppression of /? s (0). In the commen- 
surate (four-band) case, the current operator is 



7x(k) 



v x (k)T 
v x (k + Q)T 



(12) 



not. In the dSC+CSDW phase, p s (T) is exponentially 
activated: Sp s (T) oc exp(— T/S), where S is the single- 
particle gap shown in Fig. 2J In the dSC+DSDW phase, 
p s {T) is linear in T, but also has an exponential contri- 
bution from the resonance at ujq, also shown in Fig.UJ In 
most cases we have studied the exponential contribution 
is dominant. 

The inhomogeneously-doped calculations are qualita- 
tively different. The electronic potential produced by 
doping is itself weakly scattering and has little direct ef- 
fect on nodal quasiparticles 39 .; however, it disorders the 
SDW, and indirectly does have a significant effect on the 
low-w DOS. As seen in Fig. QJ characteristic features of 
the different SDW phases are washed out and N(u>) uni- 
versally obtains a dirty-d-wave form. This is one of the 
main results of this work. Not surprisingly, p a (T) also 
obtains the dirty dSC form; as shown in Fig. [5jb) and in 
the figure inset, Sp s (T) oc T 2 for T <C T c , and is linear 
in T at larger T. In general, p s {T) varies weakly be- 
tween disorder configurations. The notable exception is 
near the superconducting phase boundary (ie. w = 0.78) 
where p s (0) depends strongly on disorder configuration, 
although Sp s (T) remains quadratic in T. In macro- 
scopic samples, this will be reflected as a sensitivity to 
both sample quality and doping. The remarkable aspect 
of Fig. [5jb) is that, even for a magnetic energy scale 
M ~ 2 A (at which point p s (0) is near zero), p s (T) has 
the appearance of a dirty but pure dSC, as seen in ex- 
periments. This constitutes our main finding. 

The mechanism by which the superfluid density is de- 
pleted is quite interesting. The diamagnetic response, 
K dm , is almost independent of both w and the magnetic 



with v x — dek/dk x and tq the Pauli matrix. This matrix 
is not diagonal in the basis of Bogoliubov quasiparticles 
(ie. "f x (k) is not diagonal), meaning that the current is 
not conserved. Physically, this is because the Bogoli- 
ubov quasiparticles are formed from mixtures of states 
with crystal momenta k and k + Q. Then, because 
7a; (k) has nonzero off-diagonal matrix elements, there is a 
nonzero interband contribution in Eq. which reduces 
the overall superfluid density. We stress that this mech- 
anism is distinct from either pair-breaking or quasipar- 
ticle renormalization (ie. strong-correlation) mechanisms 
for reducing the superfluid density. 

We finish with a comment on the relationship between 
T c and p s (0). An estimate of T c as the temperature at 
which a straight line, fitted to the region T > 0.1 in 
Fig. [Ub), crosses the T-axis yields T c w 0.36, for all w. 
This is surprising, as it indicates that even for strongly 
inhomogeneous cases T c is determined by the maximum, 
rather than average, A. It also suggests that two physi- 
cal processes neglected in our calculations, phase fluctu- 
ations and glassy SDW dynamics, may play an impor- 
tant role at higher T . In particular, phase fluctuations 
are expected to be pronounced at small w where the su- 
perconductivity is spatially inhomogeneous. Glassy spin 
dynamics, provided they remain slow on electronic time 
scales, behave as quenched disordering of the SDW, and 
should not change our results qualitatively. An inter- 
esting question, outside the scope of this work, is how 
the SDW dynamics affect p s (T) at higher doping, where 
a gap in the spin-wave spectrum begins to appear. In 
summary, it seems likely that, as suggested in Refi£, T c 
is ultimately determined by a combination of phase and 
quasiparticle excitations. 
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IV. CONCLUSIONS 

In conclusion, we have shown that incommensurate 
magnetic correlations which nest the nodal points of the 
Fermi surface may coexist with <i-wave superconductivity 
with essentially no suppression of the superconducting or- 
der. Furthermore, the formation of quasistatic moments 
is a plausible explanation for the rapid suppression of 
superfluid density near p c in YBCO. We find that, pro- 
vided the spin density waves are disordered, both the 
single-particle spectrum and p s (T) are indistinguishable 



from the dirty ci-wave case. 
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